Pulling hairpinned polynucleotide chains: Does base-pair stacking interaction matter? 
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Force-induced structural transitions both in relatively random and in designed single-stranded 
DNA (ssDNA) chains are studied theoretically. At high salt conditions, ssDNA forms compacted 
hairpin patterns stabilized by base-pairing and base-pair stacking interactions, and a threshold 
external force is needed to pull the hairpinned structure into a random coiled one. The base-pair 
stacking interaction in the ssDNA chain makes this hairpin-coil conversion a discontinuous (first- 
order) phase transition process characterized by a force plateau in the force-extension curve, while 
lowering this potential below some critical level turns this transition into continuous (second-order) 
type, no matter how strong the base-pairing interaction is. The phase diagram (including hairpin-I, 
-II, and random coil) is discussed as a function of stacking potential and external force. These results 
are in quantitative agreement with recent experimental observations of different ssDNA sequences, 
and they reveal the necessity to consider the base-pair stacking interactions in order to understand 
the structural formation of RNA, a polymer designed by nature itself. The theoretical method used 
may be extended to study the long-range interaction along double-stranded DNA caused by the 
topological constraint of fixed linking number. 
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I. INTRODUCTION 

Singlc-macromolecular manipulation techniques, such 
as atomic force microscopy and optical tweezer meth- 
ods, were recently used by many authors to investigate 
the mechanical properties of double-stranded DNA (ds- 
DNA) and to understand how the tension in dsDNA in- 
fluences the interactions between DNA and proteins such 
as DNA/RNA polymerases, type II topoisomerases, and 
RecA proteins (for a brief review see Ref. [1]). Many use- 
ful insights have been obtained through a detailed anal- 
ysis of the experimental data generated by these precise 
measurements (see, for example, Refs. [2-9]). These new 
experimental techniques were recently used also on RNA 
and single-stranded DNA (ssDNA) to explore the struc- 
ture formation of these polynucleotides .BliHlij A RNA or 
ssDNA is a linear chain of nucleotide bases. At physiolog- 
ical conditions (concentration ~ 0.1M Na + and tempera- 
ture ~ 300K) , RNA molecules in biological cells can fold 
into stable native configurations (as in the case pf pro- 
teins) to fulfill their various biological functions.E3 The 
higher-order structures in RNA and ssDNA are caused 
by the possibilities to form base-pairs between the com- 
plementary nucleotide bases, A and U (T, in the case of 
ssDNA), and G and C. These structures can be further 
stabilized by the vertical stacking interactions between 
nucleotide base-pairs. In a long polynucleotide chain, 
because there are a great many different ways to form 
base-pairs, the accessible configurational space is very 
large, and the equilibrium configurations are determined 
by the competition between entropy and the interaction 
energy. To investigate how base-pairing and base-pair 
stacking influence structure formations of polynucleotide 



is therefore of both biological and theoretical interest. 

Recent measurements demonstrated that the struc- 
tures of polynucleotideSpace sensitive to the ionic strength 
of the aqueous solution.BEIrE3 In a high salt environment 
(for example, 150 mM Na + or 5 mM Mg 2+ ), electro- 
static repulsive forces between the negatively-charged nu- 
cleotide phosphate groups are largely shielded and base- 
pairing is favored; furthermore, experimental observa- 
tions suggested ihe formation of hairpinned configura- 
tions in ssDNAB To understand the hairpin formation 
in ssDNA and RNA, Montanari and Mezard suggested 
very recently a model for these polynucleotides O Base- 
pairing interaction is considered in their model, and their 
theoretical calculation can reproduce the experimental 
force-extension curve measured on a ssDNA chain whose 
nucleotide bases are arranged in relatively random order. 
It can be shown (see below and also Ref. [14]) that the 
force- induced structural transition in ssDNA is a second- 
order continuous phase transition process, characterized 
by a gradual decrease in the number of base-pairs as 
the external force is increased beyond a certain threshold 
value. 

On the other hand, Rief et a/El observed a quite differ- 
ent phenomenon: when they pulled a ssDNA chain made 
of poly(dA-dT) or poly(dG-dC) sequence apart with an 
atomic force microscope, they found that the distance 
between its two ends suddenly changes from nearly zero 
to a value comparable to its total contour length during 
a very narrow force range. 

A question thus arises: Why is the force-induced 
structural transition of ssDNA gradual and continuous 
in one case and highly cooperative in the other case, 
even though the solution conditions in these two kinds 
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of experiments are comparable? We noticed that there 
may be strong base-pair stacking interactions in the de- 
signed ssDNA chains, while stacking interactions are ab- 
sent or negligible in the relatively random sequences. 
We have previously shown that the short-ranged base- 
pair stackins-potential in dsDNA dominates the energy 
contribution^ and it is the reason both for the stability of 
dsDNA and its high extensibility at large forces Ucl Is the 
dramatically different experimental observation of Rief et 
al.t5 (compared with those of [1,11,12]) caused by base- 
pair stacking interaction? Given the fact that the stack- 
ing potential is a kind of short-ranged (nearest-neighbor) 
interaction, does the cooperative transition reported in 
[15] for such an ID system really differ qualitatively from 
that observed in [1,11]? 

To address these questions, we extend the work of 
Montanari and Mezardllil to include the effect of stack- 
ing interactions, and thus obtain a more realistic model 
for polynucleotides. We find that when the base-pair 
stacking potential is small, the ssDNA chain is in a 
hairpinncd macro-state with the base-pairs only weakly 
stacked (hairpin-I), while when the stacking potential is 
large, the ssDNA is in a hairpinned macro-state with 
almost complete pairing and stacking (hairpin-II). The 
force-induced transition between hairpin-II and random 
coil is first-order, while the transition between hairpin-I 
and random coil is continuous. There is also a continuous 
phase transition between hairpin-I and hairpin-II, which 
may be induced by changing the nucleotide arrangement 
or by environmental variations in temperature or ionic 
strength. A simplified treatment of the same question 
was reported earlier £3 

The paper is organized as following: Section [0| in- 
troduces the theoretical model of ssDNA; and Sec. |l| 
presents some qualitative arguments concerned with the 
structural transitions in the model of ssDNA; in Sec. IV 
we perform a direct comparison with experimental obser- 
vations. The main points of the paper are summarized 
in Sec. [y], accompanied by some qualitative discussions 



on RNA secondary structure formation. We suggest the 
possibility to extend the theoretical method in this pa- 
per to address the long-range interactions along topolog- 
ically constrained double-stranded DNA molecules, i.e., 
dsDNA with fixed linking number. 



II. MODEL OF SSDNA WITH BASE-PAIR 
STACKING INTERACTION 

The model used in the present study is similar with 
that developed earlier by Montanari and MezardJU albeit 
with base-pair stacking potential included. We discuss a 
fictitious polymer chain made of N tiny beads (bases) 
with index i from 1 to N (each base in the model may 
correspond to several bases in a realistic polynucleotide). 
Between any two adjacent bases (i and i + 1) there is an 
elastic bond of equilibrium length b and length variance £ 
(£ <C b, i.e., the bond is considerably rigid). For any two 
bases i and j, if their mutual distance |ry| is less than 
an interaction radius a there could be a pairing potential 
^pair( r *i) ( an d we say that an i-j base- pair is formed); 
and if any two base-pairs are nearest neighbors to each 
other (that is, i-j and (i + l)-(j — 1)), there is an ad- 
ditional stacking interaction Jf'/Zi ( r yi r i+i,j-i)- In this 
work, as wkat-is assumed in other previous theoretical 
treatments, t3~E3 only the secondary pairing patterns are 
considered, which leads to the following constraints on 
the pairing between bases: (1) each base can be unpaired 
or be paired to at most one another base; (2) if base i is 
paired to base j (suppose i < j) and k to I (k < I), then 
either i < j < k < I ("independent") or i < k < I < j 
("nested"). 

Because of these restrictions,, the partition function can 
be calculated iteratively.oEao In the case when the pair- 
ing interaction radius a is far less than bond equilibrium 
length ft (s < d), the total statistical weight Zj }i (r) for 
a polynucleotide segment (from base i to j) whose ends 
being separated by a distance r is expressed as 
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Here, fi(r) cx exp(— (|r| — b) 2 /2£ 2 ) is the prob- 
ability density for the bond vector r; /ji(r) = 
cxp[— PVp^i r ( r ij)] ~ 1 characterizes the pairing interac- 



tion between bases i and j (this viral expansion form / 
is introduced because that, in the first integral on the 
right-hand side of the equality in Eq. (oh, an vanishing 
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pairing potential between the end base j and any another 
base is assumed and the contributions of all these pos- 
sible base-paired configurations are included. Such spu- 
rious contributions should be removed since the pairing 
potential is actually nonzero); s(x) is a signal function, 
s(x) = 1 if x > and otherwise; S^f+i (Ti,j, rj + ij_i) = 



exp[— PJl'fli (rjj,rj + i >:! -_i)] — 1 characterizes the base- 
pair stacking interaction (the reason for the introduction 
of the viral coefficient g can be similarly understood as 
that for the introduction of the viral coefficient /). zj^ 
is the statistical weight for a ssDNA segment whose two 
end bases (i and j) forming a base pair: 



= s(a- |r|)exp(-/3y p J a ; r (r)) / du 1 du 2 n(u 1 )n(u 2 )Z j _ hi+1 (r - ^ - u 2 ) 
+ s ( a ~ |r|)exp(-/3V^? r (r)) J duidu 2 M(ui)M(u 2 )4i+i ( r > r - u i ~ u a) x 



Zg lji+1 (r-in-u 2 ). 



(2) 



The statistical property of this model system is deter- 
mined by the coupled equations (|l|) and (g). In the gen- 
eral case where the pairing and stacking potentials are 
dependent on base index, it is certainly of no hope to in- 
vestigate the property of the system analytically. Here we 
simplify our task by assuming that the pairing and stack- 
ing potentials between any bases have the same fcpp, i.e., 
assumingthe polymer chain to be homogeneous] ill Then 
in Eqs. (||) and (g) all the subscripts specifying specific 
bases can be dropped, especially we can write fji{v) as 



/(r), Z jti (r) as Z^r) and zjj(r) as Zf^r). In the 
model, the stacking potential is a function of the inter- 
base distances for the two base-pairs concerned. Since 
the base-pair's inter-base distance is usually very small 
(less than a), we can approximate the stacking potential 
to be constant in this range, and hence we just denote 

s(C,p) = [i + C<Kp)s(C.p)]x 

C 2 6 3 



the stacking potential to be a constant J and denote 

9 = 9ifi+l( r i,j: r i+hj-l) = const - 

The Fourier transforms of the generating functions 

(Laplace transforms) of the statistical weights are defined 

as 

S(C,P) =Jdr (jtznWCj exp(*P ■ r), (3) 



and 



S (P) (C,P) = Jdv [f^Z^^A ex P (zp • r). (4) 



Considering the iterative expressions for Z n and Zj^^ in 
Eqs. (|l|) and (g), we can show that 
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J d q( r 2 (qh(p-q)S(C,q)+.9 j dqa 2 (q) 7 '(p - q)5<"> (C, q) 



and 



aW(c, P ) = 



C 2 b 3 

(2tt) 3 



dq(T 2 (q)7'(p-q) S(C, p) + 9^ {p) (C, q) 



(5) 



(6) 



where c(p) = / dr/i(r) exp(ip • r) = 
(sinp6/p6)exp(-p 2 Z 2 /2), with <r(0) = 1; 7(p) = 
& ~ 3 i| r |<a rfr /( r )exp(ip ■ r) and Hp) = j{p) + (a/6) 3 . 
Since a -C fo, we have 7' = 7 and both 7 and 7' can be 
regarded as independent of momentum p. Then Eqs. (j^) 
and (||) lead to 



s(C,p) 



l-a(p)D(ty 



where 



D(c) = C + CS (P) (C) = -m C 3 + V2C 2 + C, 

with coefficient 771 a constant and 772 related to D: 



(7) 



(8) 



m = gl b 3 (27r)- 3 J d q( x 2 (q) = 37 (4^)- 3/2 (^), 



772(D) = D Vl 



1 + (b 2 £/gir 3/2 ) / dq 



a 2 (q) 



l l-<r(q)£> 



(9) 



(10) 



Equations (0) and @ are the central equations of this 
article. 

When a external force F is applied at the end of 
the ssDNA chain, the total partition function is Z^ = 
J drZ]y(r)exp(—(3F ■ r). The Laplace transform of this 
partition function is calculated to be 



(ii) 



N=0 
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where S(£, — i/3F) is determined by Eq. (R). 

For a linear polymer system, the free energy is an ex- 
tensive quantity proportional to the number of monomers 
N in the thermodynamic limit. This indicates that 
the smallest positive singularity point of the function 
E(£, — if3F) in the variable C corresponds to the free en- 
ergy density of the ssDNA chain.EHlcj This point is used 
in the next section. 



III. QUALITATIVE ANALYSIS OF THE 
HAIRPIN-COIL TRANSITION 

It is evident that H(£, — i(3¥) has a pole Cpoic deter- 
mined by 

£>(C P oie) - l/<7(-i/3P) - sin ^ b Fb) eM~P 2 F 2 i 2 /2), 

(12) 

where F = |F|. The function at the right side of the 
equality monotonously decreases with F from 1 at F = 
to as F — > co. On the other hand, the function D(() is 
related to the Laplace transform of the hairpinned con- 
figurations as indicated by Eq. (|J), therefore D(() has 
a finite convergence radius, and it is an increasing func- 
tion of C, before this radius is reached. The singularity of 
D(() is related to the roots of Eq. (||). Although analyt- 
ical expressions for the roots of this third-order equation 
are available, they are lengthy and here we discuss the 
statistical property of the system through an alternative 
routine. First we consider two extreme cases: 

Case A: the base-pair stacking potential J = 0. 
In this case g — and rji — and Eq. ^ reduces to 
second-order. ThiS|-situation has been studied by Mon- 
tanari and MezarcO and they found that D(() has a 
second-order branching point at Cbp equaling to the max- 
imum of the expression (—1 + \/l + 4Z3r/2)/2?72, which is 
reached at D = Db P < 1. When the external force is less 
than the threshold value F CI determined by Eq. (|l2|) at 
D = -Dbp) the polymer resides in the hairpinned phase 
with zero extension, and the free energy density equals 
to (f>(F) = (l//3)ln(bp and is force-independent. At 
F = F CT there is a second-order continuous hairpin-coil 
phase-transition (because d(/dD = at D^ p ), and the 
free energy density is changed to 4>(F) = (1//3) ln(£ po i e ). 

Case B: the stacking potential J so large that j » 1. 
In this case, Eq. ( |To| ) indicates that r]2 = Dr/i and Eq. (||) 
is equivalent to 

(C + W^)(C-£>)(C-iA/m) = o. (13) 

We readily see that when D < l/^/rji the smallest 
positive root of this equation is ( = D; and when 
D > 1/^/rji, the smallest positive solution is a con- 
stant £ = l/y/rji oc (,97) -1 ^ 2 . Most importantly, we 
have d(/dD = 1 as D approaches 1/^/fji from be- 
low. As a consequence, for F < F cr which is deter- 
mined by Eq. ( |12j ) with D = 1/y/rji, the polymer is 



in the hairpinned state and the free energy density is 
cj)(F) = -(1/2/3) In 771 oc -(1/2/3) ln( 57 ) (again indepen- 
dent of F) . For F > F cr the system is in the random coil 
state, and <j){F) = (1//3) ln(£ po i ). Here Cpoic is deter- 
mined by Eq. (|l^) (with D — Cpoic)- At F = F cr there is 
a first-order hairpin-coil phase-transition, resulted from 
the fact that d(/dD\ D= i/^- = 1. 

Comparing case A and case B, we have the impres- 
sion that the inclusion of base-pair stacking interaction 
may dramatically change the statistical property of the 
ssDNA system, even the order of the hairpin-coil phase 
transition. This is understandable qualitatively. The 
stacking interaction has two effects: (1) it makes base- 
pairing even more favorable; and (2) it causes the formed 
base-pairs to aggregate into large stacked blocks. Since 
the order of the hairpin-coil transition is related to the in- 
tensity of the stacking potential, the hairpin macro-state 
at low stacking intensity is anticipated to be different to 
the hairpin macro-state at high stacking intensity, and 
a continuous phase transition between these two hairpin 
macro-states can be predicted as the average base-pair 
stacking intensity changes. This insight is confirmed by 
observing how the order parameter of this system changes 
with stacking potential, as shown in Fig. [l]. Later we will 
refer to the hairpin macro-state at low stacking inten- 
sity as hairpin-I and that at high stacking intensity as 
hairpin-II. 

In the general case, the average number of base-pairs 
(in units of A/2) is calculated according to N^ p = 
—2d(j>/d In 7, and the average number of stacked base- 
pairs (also in units of N/2) is A s b p = —2dcj)/d\ng, 
and the relative extension of each ssDNA bond along 
the direction of the external force is obtained by Ex — 
-d<t>/dF. 

We demonstrate in Fig. |l| how the values of A^p, -A s b p , 
and the ratio A s b p /Ab p change with strength of base- 
pair stacking potential. In this figure 7 is fixed to 1.9 
(just serves as an example) and the external force is 
set to zero. We see that as the stacking potential J 
increases, all these three quantities increases and they 
approaches 1 as J reaches about J cr ~ 6ksT. The ra- 
tio iV s bp/-/Vbp could serve as an order parameter. Fig- 
ure [l] shows that when stacking potential is high, almost 
all the nucleotide bases are paired and stacked; lowering 
stacking potential from J > J cr to J < J cr there is a 
continuous phase transition where such a highly stacked 
configuration is gradually melted out and a growing frac- 
tion of bases becomes unpaired and unstacked. This is 
a structural transition between a highly stacked hairpin 
macro-state (hairpin-II) and a loosely stacked (or irreg- 
ular) hairpin macro-state (hairpin-I). This hairpin I-II 
transition is induced by changing the average effective 
base-pair stacking interaction in the polymer. Giving a 
polynucleotide chain, macroscopically the effective base- 
pairing interaction is determined mainly by its base com- 
position, while the effective stacking potential is deter- 
mined mainly by the particular arrangement of the bases 
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along the chain. Hence, for different polynucleotides with 
the same base composition, the folded hairpin configura- 
tions may be quite different and fall into two gross cat- 
alogue, dependent on their particular base sequence ar- 
rangement. This is consistent with the observation that 
certain RNA chains have unique stable native configu- 
rations. The transition between hairpins-I and -II can 
also be induced by temperature changes or variations in 
solution ionic concentrations. 

It should be mentioned that the critical stacking po- 
tential J cr is not sensitive to the pairing potential. This 
can be seen from the fact that the hairpin I-II transition 
occurs when the second term in the square brackets of 
Eq. ([To]) approaches zero, a condition which is solely sat- 
isfied when the stacking potential (and hence g) becomes 
large. 

The force-induced hairpin-coil transition is then 
second-order or first-order depending on whether the 
hairpinned configuration is type I or type II. In Fig. |^ 
the force vs extension relationship for ssDNA at different 
pairing and stacking interaction intensities are shown. As 
is expected, at large stacking intensity, there is a broad 
force-plateau and this force-plateau disappears as stack- 
ing potential is lowered. It is striking to notice that only 
the inclusion of base-pair stacking interaction can lead to 
the appearance of a highly cooperative hairpin-coil tran- 
sition. As shown also in Fig. || (the dot-dashed curve) , if 
J = the force-extension curve is always continuous and 
gradual no matter how strong the base-pairing interac- 
tion is. 

These theoretical predictions are summarized in the 
qualitative phase-diagram depicted in Fig. g[ 



IV. QUANTITATIVE COMPARISON WITH 
EXPERIMENTAL OBSERVATIONS 

In the last section we have analyzed the main qual- 
itative predictions of the present polynucleotide elastic 
model. In this section we apply these results to attain a 
quantitative understanding of some reported experimen- 
tal findings. nro 

Bustamante et al. and othersEJEillL3 have pulled a plas- 
mid ssDNA fragment of 10.4 kilobases under different 
ionic concentrations. They found that the elastic re- 
sponse of ssDNA deviates from that of a random coil both 
at high and at low ionic conditions. It was suspected that 
in high salt environments -there is partial hairpin forma- 
tion in the ssDNA chain. ErEJ As is consistent with this 
insight, Fig. ^ demonstrates that at high salt conditions, 
the ssDNA made of such a relatively random sequence 
can be well modeled by the present model with an effec- 
tive base-pairing interaction characterized by the dimen- 
sionless quantity 7=1.2 and Kuhn length b — 1.7nm and 
bond length variation £ ~ 0.0656 ~ 1.105 A. The effective 
base-pair stacking interaction for radom sequences seems 
to be quite small and is set to J = in the fitting (the 



thick solid line in Fig. [l|). Such a comparison has akeady 
been performed by Montanari and Mezard earlier .L3 For 
random sequences the reason that the effective stacking 
interaction is negligible may be explained as follows: In 
the polynucleotide chain whose bases arrange randomly, 
the formed base-pairs are usually seperated from each 
other and the possibility for two base-pairs to be adja- 
cent to each other is quite small. 

Even when the effective stacking potential J between 
two consective node-pairs becomes comparable with ther- 
mal energy fc^T (the thin doshed line in Fig. ||), the 
gradual force-extension pattern is only slightly changed. 
It indicates that the stacking interaction still plays a less 
important role than the pairing interaction. But the fit- 
ting deteriarates much for J > 2ksT if the pairing po- 
tential is kept constant (thin dotted line): for example, 
at J — 6KbT (thin solid line) there is a wide plateau in 
the theoretical force-extension curve. 

In another experiment, Rief et aL0 synthesized single- 
stranded poly(dG-dC) and poly(dA-dT) DNA chains and 
investigated their elastic responses under external force 
field. They observed a highly cooperative elongation pro- 
cess in both these two kinds of sequences, with the tran- 
sition force being 20 pN for poly(dG-dC) and 9 pN for 
poly(dA-dT). Figure |5]A and Fig. ||B shows the experi- 
mental records as well as the theoretical fittings based on 
the present model. 

In the case of poly(dG-dC), the theoretical curve is 
obtained by calculating the total extension of a polymer 
connected by a dsDNA segment of 290. Onm (with persis- 
tence length 53. Onm and stretch modulus; S = 1000. OpN 
as determined by previous experimental!) and a ssDNA 
segment of 230. Onm (with Kuhn length 1.7nm and bond 
length variance 1.105 A as determined by the data of 
Fig. ||). (We have included a segment of dsDNA simply 
because the experiment in [15] was performed by insert- 
ing a ssDNA segment between two dsDNA segments). 
In the hairpinned state at zero force, the average free en- 
ergy per Kuhn length is thus determined by the transition 
force to be <j) = -(1/2/3) 0(771) = ~5.59k B T. Similarly, 
in the case of poly(dA-dT), the fitting is obtained with 
a dsDNA segment of 80. Onm and a ssDNA segment of 
445. Onm (b = 1.7nm, I — 1.53A) with the average free 
energy per Kuhn length being <f> — — 1.763fcsT. 

Assuming each Kuhn length contains three nucleotide 
bases (i.e., each base in the model corresponds to three 
realistic nucleotides), we can infer that the free energy 
per base-pair in the poly(dG-dC) chain is 3.72kgT, while 
in the poly(dA-dT) chain is 1.18/cbjCj These values are 
close to what we estimated earlieniZl and are close to 
the phenomenological parameters Eq^c ~ 3.0teT and 
Ea-t — 1.3fcsT chosen by Bochelmann et alEB to in- 
terpret their experimental results of separating the two 
strands of a dsDNA apart by mechanical force. 



5 



V. CONCLUSIVE REMARKS 



have been observed 



s 



In this paper, we have presented an elastic model for 
single-stranded DNA and RNA polymers, where both the 
base-pairing and base-pair stacking interactions where 
considered. The theoretical calculations demonstrated 
that, depending on the intensity of base-pair stacking po- 
tential, ssDNA or RNA can form two kinds of hairpinncd 
structures, hairpin-I and hairpin-II, with the base-pairs 
in hairpin-I only slightly stacked and those in hairpin-II 
almost completely stacked. The force-induced hairpin- 
coil transition is a second-order process if the polymer 
was originally in hairpin-I macro-state and a first-order 
process if it was in hairpin-II macro-state. The phase 
diagram and the force-extension relationship for this 
polymer system were obtained and the theoretical re- 
sults achieved good agreement with experimental obser- 
vations. This work indicates the significance of base- 
pair stacking interaction to the structural property of 
polynucleotide chains. The existence of such a hairpin- 
II macro-state in a designed polynucleotide also justified 
our previous phenomenological work on the same system 
by considering only single-looped hairpin structures.E3 

The present work has the following implications. Given 
a polynucleotide chain, if the nucleotide bases are ran- 
domly arranged, then the possibility to forming stacked 
patterns of base-pairs is very small, since even forming 
the smallest stack of two base-pairs requires the correla- 
tions among four nucleotide bases. Therefore, in random 
polynucleotide chains, stacking potential is not impor- 
tant. The chains can only form loose hairpins (hairpin-I), 
and structural fluctuations between these hairpins may 
be large and quick, making it difficult for the existence of 
any well-defined stable native configurations. However, 
things may be dramatically changed if the polynucleotide 
chain composed of the same bases as the random chain 
rearranges these nucleotide bases carefully. If the bases 
in the chain is arranged in such a way that forming of 
base-pairs results in stacking of base-pairs, then it is very 
likely that the polymer will fold into certain stable na- 
tive structures which are much lower in structural energy 
than other configurations. Is this part of the reason why 
some polynucleotide sequences are folded ua^a particular 
manner (such as transfer RNA or t-RNAE3)? The nu- 
cleotide sequences in RNA molecules of biological cells 
are the results of millions of years of natural evolution 
and selection, and therefore in some sense, they should 
all be well designed. 

Finally, we just mention that the theoretical method 
used in this work to count for the long-range interac- 
tions in ssDNA may be applied to investigate the tertiary 
structures (supercoils) in double-stranded DNA caused 
by the topological constraint of fixed linking numbers. 
The topological constraint leads to a long-range inter- 
action along the dsDNA chain and make the_polymer 
to fold into compacted plectonemic structures E3 Similar 
force-extension profiles as the theoretical curves in Fig. 
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FIG. 1. The relationship between the order parameter 
Nsbp/Nbp and the effective stacking potential. Here Ab P is 
the total number of paired bases (in units of N/2), and -/V s b p 
is the total number of stacked base-pairs (also in units of N/2). 
In the numerical calculation of this figure, 7 (it accounts for 
the pairing intensity) is set to 1.9 and external force F — 0. 
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FIG. 2. The relationship between the average end-to-end 
extension of a polynucleotide and the applied force at different 
effective base-pairing and base-pair stacking intensities. Here 
b denotes the Kuhn length of a ssDNA segment, 7 accounts 
for the intensity of the pairing potential between nucleotide 
bases, and J is the effective stacking potential between nu- 
cleotide base-pairs. 
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FIG. 3. The qualitative phase-diagram of a polynucleotide 
at high-salt conditions. The phases hairpin-I and hairpin-II 
are defined in the main text. F denotes the intensity of the 
external force. 
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Relative Extension (b) 
FIG. 4. Force vs extension curves for a ssDNA polymer 
with its bases randomly ordered. The experimental data come 
from Ref. [1] (pluses) and [12] (squares), and the ionic con- 
centration is 5mM Mg 2+ . The theoretical calculation are per- 
formed with b — 1.7nm and I = 1.105A, and the other two 
fitting parameters (7 and J) are listed in the figure. 
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FIG. 5. Force vs extension curves for (A) poly(dG-dC) 
and (B) poly(dA-dT) ssDNA chains. The experimental data 
are from Ref. [15] and the theoretical fittings are done with 
b = 1.7nm, and £ = 1.1051 in (A) and £ = 1.530A in (B). The 
parameter huyi (see Eq. ^ is determined by the transition 
force, and is set to 11.18528 and 3.52741, respectively. To ac- 
count for the low extension data a segment of double-stranded 
DNA is included in the fitting (since the experiment was per- 
formed by inserting a ssDNA segment in between two dsDNA 
segments) . 
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